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Abstract 

For the local-type primordial perturbation, it is known that there is an inequality 
between the bispectrum and the trispectrum. By using the diagrammatic method, 
we develop a general formalism to systematically construct the similar inequalities 
up to any order correlation function. As an application, we explicitly derive all the 
inequalities up to six and eight-point functions. 



1 Introduction 



Primordial non-Gaussianity has been attracting attention as a powerful probe to discrim- 
inate many existing inflation models (for example, see [HE]). Statistics of zero- mean 
Gaussian fluctuations can be characterized by a variance which corresponds to two-point 
function. Hence, the non-Gaussianity of fluctuations can be linked to the higher-order 
correlation functions. A lot of studies have been done for the bispectrum, i.e. three-point 
function of the primordial perturbations. For the so-called local type perturbation, the 
bispectrum is completely characterized by a single non-linearity parameter /nl The 
current bound on / NL is — 10 < /nl < 74 at 95% confidence level [I]. 

In the last few years, trispectrum is becoming an important observable as well. For the 
local type non-Gaussianity, the trispectrum is specified by two non-linearity parameters 
tnl and #nl [HIE]- The current bounds on tnl and #nl are respectively —0.6 < tnl/10 4 < 
3.3 and —7.4 < < 8.2 at 95% confidence level [TJ. An interesting fact that boosts 

the importance of the study of the trispectrum is that tnl has a minimum determined by 



It was shown in Ref. [9] that Planck satellite can measure tnl down to 560. Therefore, 
future detection of /nl — 30 which is the central value of the current bound implies that 
we should also detect non- vanishing t N l- Detection of both / NL and r NL (and possibly <?nl 
as well) enables us to constrain or even pin down the inflation model and the origin of the 
primordial fluctuations. In fact, in Ref. [7J the authors have discussed an observational 
constraint on the ratio between / NL and tnl which is defined as A NL = tnl/(6/nl/5) 2 . 

There exist a few studies that extend the inequality (JT]). The authors of Refs. [10j ITT] 
considered one-loop corrections to flJJ) and found the one-loop effect appears as the scale 
dependent modification of the coefficient in front of f^ L . In Ref. [T2], we have discussed 
that how many parameters one needs in order to characterize higher order correlation 
functions of primordial curvature fluctuations and also introduced the non-linearity pa- 
rameters for the five-point function. The non-linearity parameters up to the six-point 
function were introduced in [T3| ITT] and two inequalities between the six- and four-point 
function were derived |13j . 

A potential importance and usefulness of the higher order correlation functions moti- 
vate us to look for the similar inequalities like ([1]) for the higher order correlation functions. 
In this paper, we provide a general formalism to construct the inequalities among the non- 
linearity parameters for the higher order correlation functions. To this end, we will adopt 
the diagrammatic approach developed in Refs. [I2J US] that turns out to be a convenient 
way to study the higher order correlation functions. By using this method, in principle, 
we can systematically derive all the inequalities up to any-point function. 

This paper is organized as follows. In the next section, we briefly review the higher 
order correlation functions of primordial curvature fluctuations, based on SN formalism 
and the diagrammatic approach. In section [21 we define the non-linearity parameters to 
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characterize the higher order functions and construct the series of inequalities by applying 
the Cauchy-Schwarz inequality. We also present all inequalities up to the six- and eight- 
point functions. Section H] is spent on the conclusion of this paper. 



2 Higher order correlation functions of 5N 

According to the SN formalism [T6l IT71 fT8l fT9] . the curvature perturbation on the uniform 
energy density hypersurface on super- horizon scales at t — tf is equal to the perturbation 
of the e-folding number evaluated from the flat hypersurface at t — t* to the uniform 
energy density hypersurface at tf at the same point: 

t(tf,x) = 6N(t f ,U;x). (2) 

The initial time t* can be chosen arbitrary. If we chose t* as the time slightly after the 
scale we are interested in leaves the Hubble horizon during inflation, then 5N in Eq. (J2J) is 
sourced by the scalar field fluctuations. Thus, £ can be Taylor-expanded in terms of the 
scalar field fluctuations as 

C(tf, X) = J2 ^N ai a 2 -aJ<p ai (U, ■ ■ ■ (U, X), (3) 

where a\, a 2 , ■ ■ ■ run from 1 to p (p is a number of light scalar fields that acquire super- 
horizon scale fluctuations during inflation). In the following, we assume that 5<p a are 
Gaussian variables. 

The two-point function of ( can be written as 

(Cfc 1 C A : 2 ) = (27r) 3 ^ 1 + fc 2 )P c (fc 1 ). (4) 

To leading order in the field fluctuations, (k) is given by 

P € (k) = N a N a P(k), P(k) = — (j£J , (5) 

where summation over a is assumed. 

The three-point function of ( can be written as 

(CfeCfeCfi.) = (2n) 3 5(h + k 2 + h)B ( (h, k 2 , h). (6) 
To leading order in the field fluctuations, is given by 

AT NkN u 

B c (h, k 2 , fe s ) = ( P dh)Pdh) + 2 perms.) . (7) 

Following literatures, the constant coefficient in front of the square of the power spectrum 
is written as [20] 

, 5 N a N b N ab 

/NL " 67W (8) 
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The four-point function of £ can be written as 



(CfcCfcCftfo) = (2*fS(h + h + k 3 + h)T c (h, k 2 , fc s , fc 4 ). (9) 
To leading order in the field fluctuations, Tq is given by 

T c (fc 1; k 2 , k 3 , k A ) = N " N "f T b f C (P c (ki)Pdki2)Pdh) + 11 perms.) 

+ N * N J N A C Z abC (Pdh)Pdk2)Pdh) + 3 perms.) , (10) 

where ky = \ki + kj\. Note that there appear two distinct terms that exhibit different 
wavenumber dependence. As a consequence, we need two parameters to specify the tri- 
spectrum. Following literatures, the two constant parameters are defined by [6] 

N a N ab N bc N c 25 N a N b N c N abc 

TNL " (W ' ^ NL ~~ 54 {NdNd f • ^ 

We can further proceed to higher order correlation functions as follows. The leading 
order of the n-point function is given by a sum of several distinct terms which are products 
of (n — 1) power spectra and exhibit different wavenumber dependence. According to the 
diagrammatic method [12] , each of these leading terms has a corresponding connected tree 
diagram that consists of n vertices and (n — 1) lines connecting two vertices. Tree diagrams 
for the two-, three- and four-point functions are shown in Fig. [U In reverse, given a tree 
diagram with n vertices, we can reconstruct the corresponding term that constitutes the 
n-point function as follows. First, we assign a different wavenumber {k\, • • • , k n } to each 
vertex of the diagram, where {k±, ■ ■ ■ ,k n } are the arguments of the n-point function with 
the constraint k\ + - ■ -+k n — 0. Next, we assign a wavenumber to each line in the diagram, 
too. Removing a line from the diagram yields two respectively connected sub-diagrams. 
Then, one assigns to the removed line the sum of the vectors associated with all vertices 
in one of the two sub-diagrams. Any use of the two sub-diagrams yields the same answer 
because of the constraint k± + • • ■ + k n = 0. 

After associating the wavenumbers with all lines, we can assign the corresponding 
factors to the vertices and the lines. As for the vertex with p lines attached, assign the 
factor N a i... a to it. As for the lines, assign P, where the argument of P is set to the length 
of the wavenumber associated with each line. By multiplying all these factors assigned 
to vertices and lines, we obtain the corresponding leading term. Any term constructed 
in this way is given by a product of two terms. One is the constant term which is given 
by contracting all the pairs of indices of the expansion coefficients assigned to each vertex 
of the tree diagram. The other term is the (n — 1) product of the power spectra of 
the scalar field plus its permutations. If we replace each power spectrum of the scalar 
fields by that of the curvature perturbation, i.e. P —> P^/(N a N a ), the constant part 
gets a factor (A^ a A^) n_1 in the denominator. Following the cases for the bispectrum and 
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Figure 1: The tree diagrams corresponding to the power spectrum, bispectrum and trispec- 
trum, based on 5N formalism. We find that we need two parameters in order to charac- 
terize the trispectrum. 

the trispectrum, we define the constant parameter corresponding to this diagram as a 
coefficient in front of the (n — 1) product of P^. By taking the sum over all the possible 
(connected) tree diagrams with n vertices, we obtain the n-point function. 

It was shown in Ref. [12] that the functions constructed from two tree diagrams that 
are not isomorphic each other always yield different wavenumber dependence. Therefore, 
the number of independent constant parameters required to specify the n-point function 
is equal to the one of all the possible connected tree diagrams with n vertices that are not 
isomorphic to each other. 

A formula for a number t n of the independent connected tree diagrams with n vertices 
is given in Refs. [2T|[22]. According to Refs. [211 [22], we can formally construct a function 
t(x) as an infinite series, 

t(x) = J2 t nX n - (12) 
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Then t(x) is given by 

I i 

t(x) = r(x) - -r 2 (x) + -r(x 2 ), (13) 
where r(x) is a function that satisfies, 

(14) 

From the last equation, we can recursively obtain the Taylor-expansion coefficients of r(x). 
Then from Eq. f|T3|) . we can recursively obtain t n as well. For example, t 2 = 1, t 3 = 1, t 4 = 
2, t 5 = 3, t 6 = 6, t 7 = ll, t 8 = 23,---. 

3 Inequalities of the Cauchy-Schwarz type 

As we mentioned in the last section, wavenumber dependence of the n-point function is 
completely specified by a set of constant parameters each of which has a corresponding 
connected tree diagram with n vertices. For the purpose of constructing the general 
inequalities among the non-linearity parameters, it is convenient to adopt a normalization 
condition, 

N a N a = 1. (15) 

Without a loss of generality, this condition can be always imposed by rescaling the scalar 
field fluctuation. With this condition, the denominator of the non-linearity parameters for 
the n-point function, which is given by (N a N a ) n ~ , becomes unity. 

3.1 Definition of the non-linearity parameters 

For any connected tree diagram with n vertices, if we cut a line appearing in the diagram, 
the diagram splits into two sub-diagrams. If the resulting two sub-diagrams are isomorphic 
to each other, we call the parent diagram symmetric diagram. Obviously, in order for a 
diagram to be a symmetric diagram, n must be even. Let us denote by K2 m a number 
of all the symmetric diagrams among t2 m tree diagrams with 2m vertices. For example, 
^2 = 1, ^4 = 1, ^6 = 2, k 8 = 4, Kio = 9, • • ■ . Our aim is to construct the inequalities 
of the Cauchy-Schwarz type for the non-linearity parameters up to the n-point functions. 
As it will be clear later, this is possible when n is even. For the case of odd n, we can 
not construct the inequalities. For example, there are inequalities among the non-linearity 
parameters up to the six-point function. But there are no inequalities in which the highest 
order of the non-linearity parameters is five. 

Since whether n is even or odd is important, let us first write the non-linearity param- 
eters for the 2m-point function as 

rpf 1 ) p(*2m)x _ J J (1) (K2mh [A 1 ) (*2m-K2mhl 

l r 2mi 5 r 2m J ~ lV2m! * * ' ) '2m h li^m' ' ' ' 1 Vim > > ' \ LV> ! 
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Figure 2: The symmetric and non-symmetric diagrams for n = 4, 6 cases. 
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Figure 3: The diagrams for n = 3, 5 cases. 

Here, r 2 ^ and g^m are the non-linearity parameters for the symmetric diagram and the 
non-symmetric diagram, respectively In Fig. [21 we show the symmetric diagrams and 
the non-symmetric diagrams for n = 4, 6-point functions. For the four-point function, a 
diagram corresponding to tnl is a symmetric diagram while the one corresponding to #nl 
is a non-symmetric diagram. Therefore, relations between (74 , g±) and the conventional 
non-linearity parameters tnl, S'nl are given by 



tnl, 



9? 



54 
25 



<7nl- 



(17) 



For (2m — l)-point function, we write the non-linearity parameters as 



r 2m J — li2m-H 



2m-l 



In Fig. El we show the diagrams for n = 3, 5-point functions. A relation between /g and 
the conventional non-linearity parameter /nl is given by 



(1) 



6 



/] 



NL- 



(19) 



3.2 Construction of the inequalities 

Let us suppose that the highest order of the non-linearity parameters appearing in the 
inequalities is n = 2m. We can construct the series of inequalities in a following way. 
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3.2.1 The highest order n = 2m 

Since is a non-linearity parameter corresponding to the symmetric diagram, it can be 
written as 



T. 



2m 



V£ ] ■ (20) 
where Vm is a vector associated with a diagram Dm which is one of the sub-diagrams 



generated by dividing the symmetric diagram into two diagrams that are isomorphic to 
each other. Therefore, Dm is a tree diagram with m vertices and with one external line. 
Next, let us consider a quantity 

yW.f(i) (21) 

for % 7^ j. This must be a non-linearity parameter whose corresponding diagram is a 
diagram with 2m vertices and is given by connecting the external line of Dm with that of 
Dm . Since this diagram is a non-symmetric diagram, there exists some k such that 

v m l) ■ v$> = g&. (22) 

Applying the Cauchy-Schwarz inequality to Eqs. (1201) and fl22l) yields an inequality, 

r 2m r 2m — \92m) ■ (23) 



A number of the independent inequalities of this kind is K 2m (K 2m . — 1) /2. These are the 
inequalities among the non-linearity parameters for 2m-point function. 

3.2.2 Next highest order 2m — 1 

Let us then lower the order by one and consider a quantity, 

V®-V£ v (24) 

This is a non-linearity parameter whose corresponding diagram is obtained by connecting 
the external line of D m with that of D^_ x . Therefore, there exists some k such that 

Vm ' Kra-l = fim-l- (25) 

Meanwhile, since a quantity Vm-i ' K^-i ^ s a non-linearity parameter whose corresponding 
diagram is obtained by connecting the external lines of the same two D^_ x , we have 

^•^1=^. (26) 
Applying again the Cauchy-Schwarz inequality to Eqs. ( |25l) and ( )26l) yields 

(0 Ji) > f f (*) V (27) 

A number of the independent inequalities of this kind is K,2m X ^2m-2- These are the 
inequalities among the non-linearity parameters for 2m, 2m — 1 and 2m— 2-point functions. 
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3.2.3 Next-to-next highest order 2m — 2 

Let us further lower the order by one and consider a quantity 

■ V^_ 2 . (28) 

This is a non-linearity parameter whose corresponding diagram is obtained by connecting 
the external line of D m with that of D m ^_ 2 . Since this diagram can be both symmetric or 
non- symmetric, we write it as 

Applying the same reasoning as before, we obtain 



T 2m r 2(m-2) — (^Li-^) ■ (30) 

A number of the independent inequalities of this kind is K 2 m x K 2m _4. These are the 
inequalities among the non-linearity parameters for 2m, 2(m — 1) and 2(m — 2)-point 
functions. 

By repeating the above procedures until the order reduces to one, we can obtain all 
the Cauchy-Schwarz type inequalities that involve the non-linearity parameters for the 
2m-point function as the highest order. In particular, at the final procedure, we end up 
with getting a set of inequalities: 



r£>(/i fc |i) 2 - (31) 



A number of the independent inequalities of this kind is ^m- 

Combining all these results, a number a^m of independent inequalities that involve the 
non-linearity parameters for the 2m-point function as the highest order is given by 

&2m = K 2m (^{.^m ~ 1) + K 2m-2 + K 2m-4 + ' ' ' + K 2^j ■ (32) 

For example, = 1, = 5, a 8 = 22, aio = 108, 



3.3 Application to four, six and eight-point functions 
3.3.1 Case of four-point function 

There are two different tree diagrams for the four-point function. One of them is a sym- 
metric diagram which corresponds to tnl and the other one is a non-symmetric diagram 
corresponding to #nl- According to the general argument we developed in the last sub- 
section, there is no inequality between tnl and g NL because k 4 = 1. Then, let us lower the 
order by one. At this stage, we obtain one inequality 

74 > fl (33) 
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In terms of the conventional non-linearity parameters, this can be written as 

7-nl > ^fL- (34) 
Since we have reduced to the lowest order, there are no other inequalities. 

3.3.2 Case of six-point function 

Since Kq = 2, there exists a single inequality among the non-linearity parameters for the 
six-point function: 

4 1} 4 2) > (#) 2 . (35) 

where the correspondence between the diagrams and the non-linearity parameters is given 
in Fig. [2j By lowering the order by one, we obtain k 6 x k 4 = 2 inequalities: 

r<"r 4 > (/<") 2 , rf r 4 > (/f ) 2 . (36) 
Reducing further the order by one, we obtain k 6 x k 2 = 2 inequalities: 

r 6 (1) > rl rf > gj. (37) 

The last two inequalities were also provided in [T3]. Since we have reduced to the lowest 
order, there are no other inequalities. Total number of inequalities is five. 

3.3.3 Case of eight-point function 

Since k$ = 4, there exist six inequalities among the non-linearity parameters for the eight- 
point function: 



r 8 (1) rf > (rf>) , r<"rf > (#) , 



T M-C) > fjl)" 1 2 
I a ' s i_ 



rfxf > ( 9 «) 2 , tPtP >(£>)', rf r« > ( 9 f ) 2 . (38) 

where the correspondence between the diagrams and the non-linearity parameters is given 
in Fig. HI By lowering the order by one, we obtain k$ x k 6 = 8 inequalities: 

rfrf > (/f) 2 , r«"rf > (jf>)', rf r<" > (/»)" ', rf rf> > (#>) 2 , 

rfV«'> > (/f ) 2 , rfrf > (/<<*> )', xf x« > T « T f > (/<» " 

where the correspondence between the diagrams and the non-linearity parameters is given 
in Fig. [51 Reducing further the order by one, we obtain Kg x = 4 inequalities: 

r s (1 V 4 > ( S «) 2 , rfr 4 >( 9 f) 2 , rfr 4 > ( 9 f ) 2 , „ > (r«) 2 . (39) 
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Figure 4: The correspondence between the diagrams and the non-linearity parameters for 
8-point function. 
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Figure 5: The correspondence between the diagrams and the non-linearity parameters for 
7-point function. 
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Reducing further the order by one, we obtain k§ x k 2 = 4 inequalities: 




Since we have reduced to the lowest order, there are no other inequalities. Total number 
of inequalities is 22. 

4 Conclusion 

The primordial non-Gaussianity has been focused on by many authors as a new probe of 
the inflation dynamics. In addition to the bispectrum, use of the higher order correlation 
functions will become useful in the future. In particular, the inequality between tnl 
and (?nl shows the importance of the trispectrum to look for the non-Gaussianity in the 
primordial perturbations. 

In this paper, we developed a general formalism to construct the Cauchy-Schwarz type 
inequalities among the non-linearity parameters for the higher order correlation functions. 
This method enables us to derive all the inequalities up to any-point function. As an 
application, we explicitly derived all the inequalities up to the four, six and eight-point 
functions. We first confirmed that there is just one inequality, which is Eq. (JTJ), up to the 
four-point function. Up to the six-point function, there are five new inequalities two of 
which were given in Ref. [13]. Up to the eight-point function, there are further 22 new 
inequalities. 
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